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Abstract 


A  two-scale  approach  based  on  eigendefonnation-based  homogenization  is  explored  to  predict 
the  behavior  of  concrete  targets  subjected  to  impact  loading  by  high  speed  projectiles.  The 
method  allows  accounting  for  micromechanical  features  of  concrete  at  a  computational  cost 
comparable  to  single  scale  phenomenological  models  of  concrete.  The  inelastic  behavior  of 
concrete  is  modeled  using  three  types  of  eigenstrains.  The  eigenstrains  in  the  mortar  phase 
include  pore  compaction  (or  lock-in),  rate-dependent  damage  and  plasticity  eigenstrains,  whereas 
the  inelastic  behavior  of  aggregates  is  assumed  to  be  governed  by  plasticity  only.  Material 
parameters  were  identified  using  inverse  methods  against  unconfined  compression  and  uniaxial 
compression  tests.  A  unit  cell  was  constructed  from  a  3D  digital  image  of  concrete.  The 
eigendeformation-based  homogenization  approach  was  validated  for  the  projectile  penetration 
into  concrete  target.  The  simulation  results  were  found  to  be  in  reasonable  agreement  with  the 
experimental  data.  Attention  is  restricted  to  non-reinforced  concrete. 

1.  Introduction 


Portland  concrete  has  been  widely  used  as  a  construction  material  for  civil  and  military 
applications  since  the  18th  century,  even  though  its  grandfather,  Roman  concrete  made  from 
quicklime,  pozzolanic  and  an  aggregate  of  pumice,  dates  back  to  Roman  Empire.  The  reasons  for 
concrete’s  ubiquity  in  practice  are  not  only  it  being  one  of  the  environmentally  friendliest  and 
cost-effective  materials,  but  also  due  to  its  high  compressive  strength,  which  makes  it  play  a 
critical  role  in  defense  applications  including  nuclear  reactor  containments  and  fortification 
installations. 

The  problem  of  predicting  the  behavior  of  concrete  targets  subjected  to  impact  loading  by 
high  speed  projectiles  has  been  of  military  interest  for  many  years.  Various  procedures  have  been 
developed  and  they  generally  fall  into  three  categories:  (1)  empirical  or  analytical  approaches 
[1,2],  (2)  experimental  approaches  [3-5]  and  (3)  numerical  approaches  [6-15]. 

For  numerical  approaches,  such  as  the  finite  element  method,  to  be  viable,  it  is  essential  to 
characterize  the  constitutive  models  for  concrete.  Concrete  is  typically  modeled  as  a 
homogeneous  single-scale  [7-11,16,17]  or  heterogeneous  two-phase  [6,12-14]  material.  In  an 
attempt  to  capture  the  dilatant  behavior  of  concrete  the  cap  plasticity  model  was  employed  in 
[7,8],  In  [9],  a  modified  Hohnquist-Johnson-Cook  (HJC)  model  for  concrete  with  enhanced 
pressure-shear  behavior,  strain-rate  sensitivity  and  additional  damage  variables  tracking  the 
tensile/shear  cracking  and  pore  compaction  were  introduced.  A  coupled  plasticity-damage  model 
was  considered  in  [10],  whereas  an  energy  release  rate  based  plasticity-damage  model  was 
suggested  in  [1 1].  In  [17]  a  microplane  model  of  concrete  was  developed  to  account  for  pressure 
sensitivity,  dilatancy,  deviation  from  nonnality,  Bauschinger  effect  and  hysteresis. 

Consideration  of  concrete  as  a  heterogeneous  or  composite  material  consisting  of  mortar 
and  aggregates  has  been  somewhat  less  popular  (see  [6,12-14]  for  noteworthy  exceptions),  The 
premise  of  this  approach  is  that  each  phase  can  be  modeled  using  different  material  laws  and  the 
interaction  between  the  phases  can  be  explicitly  accounted  for.  This  is  of  particular  importance 
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since  elastic  degradation  of  mortar  exhibits  a  more  compounded  process  than  that  of  the 
aggregate.  In  the  mortar  phase,  microcrack  closure  under  reversal  loading  gives  rise  to  stiffening 
of  material  and  microcrack  growth  is  not  perfectly  brittle  but  rather  involves  plastic 
deformations.  On  the  other  hand,  the  aggregate  phase  can  be  adequately  described  within  the 
framework  of  classical  plasticity.  The  interaction  between  the  phases  so  far  have  been 
predominantly  accounted  for  by  means  of  various  effective  medium  models,  such  as 
Mori-Tanaka  and  self-consistent  approaches.  Notable  exceptions  were  reported  in  [18-20]  where 
computational  homogenization  methods  were  employed  to  explicitly  account  for  concrete 
microstructure.  Also  noteworthy  are  various  mesomechanical  models  where  continuum  is 
replaced  by  discrete  medium  [21,  22], 

While  computational  cost  of  the  computational  homogenization  approaches  is  a  small 
fraction  compared  to  the  direct  numerical  simulation,  where  a  characteristic  mesh  size  is  of  the 
heterogeneity  order,  they  remain  computationally  prohibitive  for  complex  microstructures  such 
as  those  of  concrete.  This  is  because  a  nonlinear  unit  cell  problem  for  a  two-scale  problem  has  to 
be  solved  for  a  number  of  times  equal  to  the  product  of  the  number  of  quadrature  points  at  a 
macroscale  and  the  number  of  load  increments  and  iterations  at  the  macroscale.  If  this  approach 
is  to  be  a  viable  alternative  to  modeling  concrete,  it  is  necessary  to  reduce  the  computational 
complexity  of  the  fine  scale  (unit  cell)  problem  without  significantly  compromising  on  the 
solution  accuracy  of  interest,  which  is  typically  at  the  macro  level.  In  this  manuscript  we  explore 
the  application  of  the  eigendefonnation  based  homogenization  [23-28]  to  multiscale  modeling  of 
the  non-reinforced  concrete  targets  subjected  to  impact  loading  by  high  speed  projectiles.  The 
salient  feature  of  the  eigendeformation-based  homogenization  is  that  the  unit  cell  problem  is 
formulated  in  terms  of  eigendefonnation  modes,  which  a  priori  satisfy  equilibrium  equations  in 
the  microscale,  and  thus  eliminate  the  need  for  costly  solution  of  discretized  nonlinear 
equilibrium. 

The  paper  is  organized  as  follows.  The  problem  statement  is  formulated  in  Section  2.  The 
formulation  of  eigenstrains  and  the  overall  solution  algorithm  are  given  in  Section  3.  Validation 
studies  are  considered  in  Section  4  starting  with  the  digital  image  of  the  concrete  microstructure, 
calibration  of  microphase  properties  and  prediction  in  a  penetration  problem.  Finally,  several 
concluding  remarks  are  drawn  in  Section  5.  The  eigendeformation-based  homogenization  theory 
is  outlined  in  the  Appendix. 
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2.  Problem  Statement 


Figure  1.  Macro-  and  micro-  structure  of  concrete  material 


Consider  a  concrete  as  a  heterogeneous  material  formed  by  a  repetition  of  locally  periodic 
microstructure  as  shown  in  Figure  1.  The  microstructure  consisting  of  mortar  and  aggregate 
phases  is  described  by  a  unit  cell  or  so-called  representative  volume  element  (RVE). 

The  macroscopic  problem  considered  in  the  paper  is  a  penetration  example  [4]  by  a  rigid 
projectile  as  shown  in  Figure  2.  The  projectile  geometry  is  given  in  Figure  3.  The  physical 
properties  of  the  concrete  used  in  this  paper  are  listed  as  Table  1. 

The  governing  equations  describing  the  response  of  the  structural  system  are: 

Momentum  balance: 


o-jj  (x,  t )  +  b:  (x,  t)  =  p(x,  t)ih{x,  t). 


(1) 


where  cr.  is  stress,  h  is  body  force,  p 
Kinematics  relation: 


w  ,1 
£ij(x’t)  =  - 


dui  dti; 
'  +  ] 


dxj  dx( 


is  density  and  u. 


is  displacement. 


(2) 


where  ui,sj  are  displacement  and  strain  components,  respectively.  The  superimposed  dot 

denotes  material  time  derivative. 

Constitutive  relation: 


aj(xJ)  =  Ljkl(x) 


£k /Cm)-£  V«Cm) 

i 


(3) 


where  Lijkl  is  elastic  moduli  and  1  p:j  is  the  I- th  component  of  the  eigenstrain  We  assume  an 

additive  decomposition  of  total  strains  into  elastic  strains  and  eigenstrains  consisting  of  inelastic 
strains  (damage,  plasticity,  etc.),  hygrothemal  strains,  and  phase  transformation  strains.  The 
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formulation  of  the  inelastic  deformation  strains  is  detailed  in  Section  3. 


height=l  .22m 


diameter=1.83m 

Figure  2.  Definition  of  the  penetration  model 


Dia.  80 Cl  (3. 150) 


Figure  3.  The  geometry  of  3.0  CRH  projectile  |4] 
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Table  1.  Physical  properties  of  concrete  used  in  the  paper 


Compressive  strength 

38.7MPa 

Wet  density 

2.25  Mg/ m3 

Volume  solids 

82.53% 

Volume  water 

5.63% 

Volume  air 

11.84% 

Volume  voids 

17.47% 

The  mixture  portions  of  the  concrete  (38.7MPa)  are  listed  as  Table  2. 
Table  2.  Mixture  portions  of  concrete  used  in  the  paper 


Materials 

Mass  in  pounds 

ASTM  Type  I,  Holnam  Cement  Co. 

543 

Natural,  siliceous,  concrete  sand,  SSD  condition 

1478 

ASTM  No. 89,  3/8  NMS,  crushed  limestone,  SSD  condition 

1566 

Water  reducer,  high-range,  Eucon  37 

961  ml  or  32.5  fl  oz 

Water  reducer,  nonnal-range,  Eucon  WR-91 

961  ml  or  32.5  fl  oz 

Water 

315 
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3.  Microscale  Inelastic  Properties  of  Concrete:  Eigenstrain 


The  theoretical  framework  of  the  eigendeformation-based  approach  is  outlined  in  Appendix 
A.  In  this  Section  we  focus  on  the  formulation  of  the  inelastic  deformation  eigenstrains  1  n 

assuming  perfect  interfaces,  i.e.,  neglecting  the  eigenseparations  §h .  Having  been  frequently 
referred  to  as  plastic-fracturing  materials,  concrete  exhibits  plastic  flow  in  both  mortar  and 
aggregate  phases.  We  hereby  define  plastic  strain  sp  as  one  source  of  the  eigenstrains  in  both 
phases.  Another  source  of  eigenstrains  denoted  as  sd  is  formed  by  crack  nucleation  and 
microcracks  growth,  which  typically  takes  place  in  the  mortar  phase  [13].  Finally,  the  closure  of 
microcracks  and/or  original  air  voids  under  compressive  load  gives  rise  to  the  compaction  of 
concrete  material.  This  causes  a  sudden  stiffening  of  the  mortar  phase  [9,13]  known  as  lock-in. 
This  phenomenon  in  the  mortar  phase  will  be  modeled  by  so-called  lock-in  eigenstrain  denoted 
as  sc .  The  next  two  subsections  will  focus  on  the  formulation  of  the  eigenstrains  in  the  two 
phases.  The  algorithm  details  will  be  given  in  Section  3.3. 

3.1.  Eigenstrains  in  the  aggregate  phase 


The  inelastic  behavior  of  aggregates  is  modeled  using  classical  isotropic/kinematic 
plasticity  theory.  Thus  plastic  strain  is  assumed  to  be  the  only  eigenstrain  source 

/J(agg-)  =£P{agg\  (4) 


The  inelastic  parameters  for  the  aggregate  phase  are:  hardening  modulus  #(ogg) ;  yield  stress  as 
ct^88^ ,  and  mixture  parameter  0{ag8)  defining  the  relative  weight  of  the  isotropic  and  kinematic 

hardening  (  0{ag8  )  =  1  for  pure  isotropic  hardening  and  Qu‘Kg:)  =  o  for  pure  kinematic 
hardening).  Two  additional  internal  variables  denoted  [a,/)]  are  defined  with  a  being  the 
equivalent  plastic  strain,  and  b  tracking  the  center  of  von  Mises  yield  surface  in  the  stress 
deviator  space.  Given  local  strain  in  the  aggregate  sfss:),  one  can  find  local  plastic  strain 

£p(agg-">  by  racijal  return  or  similar  algorithm  [29]. 


3.2.  Eigenstrains  in  the  mortar  phase 


We  assume  that  the  three  sources  of  eigenstrains  £p(monar}  ?  £d(mortar) ,  £c(monar)  [  ntrodUccd 
in  Section  3  may  not  evolve  simultaneously.  The  damage  eigenstrain  phmortan  is  assumed  to  be 
initiated  and  accumulated  under  tensile  loading  at  the  microscale,  while  £ pl mo"ar '  and 
£c(mortar )  evo]yg  unc[er  compressive  loading.  The  tensile  and  compression  loadings  are  determined 
by  probing  the  macroscopic  volumetric  strain  5  =  £,  ,  where  Sl  is  a  macroscopic 

principal  strain.  The  evolution  of  the  inelastic  deformation  eigenstrain  for  mortar  is  then 
expressed  as 
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^-j-^mortar)  _  ^p(mortar)  _|_  £C {mortar) 

j  jj{mortar)  _  {mortar) 


//S<0, 

ifZ>0. 


(5) 


For  plastic  strain  £p(mortar)  ?  we  employ  the  same  linear  isotropic/kinematic  hardening  law  as  for 
the  aggregate  phase.  The  formulation  of  the  other  two  eigenstrain  components  in  the  mortar 
phase  is  discussed  next. 

3.2.1.  The  damage  eigenstrain  £d<mortar'> 


We  consider  a  rate-dependent  damage  model  to  alleviate  (at  least  partially)  mesh  size 
dependency  caused  by  strain  softening.  For  simplicity,  we  omit  the  superscript  (mortar).  The 
damage  parameter  <x>  is  defined  to  be  a  function  of  equivalent  strain  i ,  i.e.  <x>  =  O(f') .  The 

damage  eigenstrain  is  defined  as  sd  =  cos .  With  the  internal  variable  r ,  which  records  the 
largest  equivalent  strain  in  loading  history,  we  define  the  closure  of  damage  function  as 

C.={(£,r)ef?xf?|g(i,r)<0}  (6) 

where  g{e,r )  :=s  —  r.  The  damage  evolution  is  defined  as 


•  1  /  t*  x\d®(£) 

<o  =  -{g{e,r))- 


& 


ds 


r  =  ^{g{£,r)) 


where  {■)  is  the  Macaulay  bracket  and 


x)  = 


(x+  I  X I) 


(7) 


.  The  rate  effect  is  controlled  by  <9 .  If 


3  — >  0 ,  the  model  becomes  rate-independent,  whereas  if  ,9  cc  ,  the  model  coincides  with  the 
instantaneous  elastic  response.  Various  functions  of  damage  parameter  <x>  =  ©(£)  have  been 
proposed.  Here,  we  adopt  the  simplest  variant  defined  in  Figure  4. 

In  Figure  4,  S  is  the  critical  stress  at  the  end  of  elastic  process,  G  is  the  area  under  the 
stress/strain  curve  or  the  strain  energy  density.  The  damage  parameter  <x>  is  defined  as 


co  =  0(£)  = 


O) 


0 


kf-k, 


a  <  kt 
k<s<k 
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£>kf 


_-S  +  ^24GEa-3S2 
"  2 Ea  ’ 

where  E  is  elastic  modulus  in  the  loading  direction. 


E 


(8) 


(9) 


7 


<i> 


Figure  4.  Relation  between  cumulative  damage  function  and  stress/strain  curve 


3.2.2.  The  lock-in  eigenstrain  a 


c  {mortar) 


The  lock-in  eigenstrain  s'  is  assumed  to  be  a  function  s'  =  KY(a)  of  the  characteristic 
compressive  strain  defined  as 

£  =  Z(-^)’  (i°) 

1=1 

where  (•)  is  the  Macaulay  bracket  and  £:  is  local  principal  strain.  The  function  T*  is 
characterized  by  three  material  parameters  p.,  pf  and  s‘max  as  (see  also  Figure  5) 

0  a  <pt 

I  £max(e  ~Pi) 


£c=V(s)  = 


Pf-Pi 


P^£  ^Pj 
£  ^ Pf 


(ID 
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Figure  5.  Relation  between  and  characteristic  strain  £ 


3.3.  Numerical  implementation 

Recalling  the  internal  variables  a ,  ^ ,  and  r ,  and  defining  the  internal  variable  s  to 
record  the  largest  characteristic  compressive  strain  £  in  the  loading  history,  the  algorithm  in 
Box  1  describes  an  iterative  scheme  for  updating  multiple  eigenstrains.  The  algorithm  takes 
advantage  of  the  reduced  order  equation  (A  10)  given  in  the  Appendix.  Boxes  1.1 -1.3  give  the 
corresponding  subroutines  for  updating  different  eigenstrains. 
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Box  1.  Main  Algorithm:  Updating  multiple  eigenstrains  in  concrete  (the  reduced  order  unit  cell  problem) 

Given  microscopic  variables  (subscripts  denoting  the  time  step): 

Microscopic  strain  £J-mo"ar'> ;  £(asg  ) 

Plastic  eigenstrain  £P(monar ) ;  £p(asg-) ,  an(j  piastic  internal  variables 
Damage  eigenstrain  £l^moru"  1  ,  ancj  damage  internal  variables 
Lock-in  eigenstrain  £^mor,a'  >  ,  an(j  Jock-in  internal  variable 
and  macroscopic  strain  s  at  time  tn  ,  strain  increment  A£n+l  between  tn  and  tn+l  . 

Find  the  above  micro-  and  macroscopic  variables  at  time  t  x . 

1 .  Let  p  ,  —  p  “F  A p  .  A?  —  t  i  —  t  . 

Cjt  1  bn  b  n+l  ’  n+ 1  n 

2.  Calculate  the  macroscale  principal  strain  £l  . 

3.  Calculate  the  macroscale  volumetric  strain  5  ±l  =V  p,  . 

n+\  Z^I=\bIll+l 

4.  Using  £Jmol‘ar)^£Jagg^  as  initial  values,  apply  Newton  Method  to  solve  for  s\'l°rtar)  from 

the  equations  (see  equation  (A9)  in  the  Appendix) 

( mortar )  p(mortar., mortar)  f  p(mortar)  .  d(mortar)  .  c(mortar)  |  -p(mortar.,agg.)  p(agg.)  Amortar) — 

V (n+l)  ykl  y  kl(n+\)  kl(„+ 1)  kl(n+\)  )  Ukl  kl(„+\)  ijkl  £  kl(n+ 1) 

c-(agg-)  _  p(agg.,agg.)  p(agg.)  _  p(agg mortar)  (  p(mortar)  .  d(mortar)  .  c(mortar)  |  _  Aagg)~ 

U  (n+l)  ykl  kl(n+ 1)  ijkl  y  kl(n+ 1)  kl(„+ 1)  kl(n+ 1)  J  ijkl  £  kl(n+ 1) 

For  each  iteration  of  the  Newton  Method: 

(a)  update  the  eigenstrains  £P^morm  >  ancj  £^Jas"  ]  using  the  subroutine  in  Box  1.1. 

(b)  update  the  eigenstrain  £^Jmo,,,n  1  using  the  subroutine  in  Box  1.2. 

(c)  update  the  eigenstrain  £(nJm°"‘"  1  using  the  subroutine  in  Box  1.3. 

5-  Let  (•)„  =  (')«+!• 
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Box  1.1.  Subroutine  for  updating  plastic  eigenstrain  £p+l 

Given  eigenstrain  £p  ,  plastic  internal  variables  an  ,  bn  ,  microscopic  strain  £n+l  ,  and  macroscopic 
volumetric  strain  S;|+1 . 

Find  £n+l ,  an+l ,  bn+l . 

If  Sn+1  >  0,  then  (tensile  loading): 

<+!=<’  a„+l=an’  K+ l=bn- 

Else  (compressive  loading): 


(a)  find  the  deviator  en+1  =  £n+l  -  -  ( tr[sn+l  ])/  . 

(b)  compute  the  trial  stress  =  2 jU(en+l  -  ep  )  -  bn  . 

(c)  check  the  yield  condition  fn+]  :=||  Cnlf  II  — V 2/3  (<JY  +OHan  )  .  if  fn+x  >  0 ,  go  to  next 


else  A/n+l  =  0 ,  go  to  step  (f). 

(d)  calculate  the  consistency  parameter  A yn+1  = 


2jU  +  ^H 


(e)  calculate  the  unit  vector  normal  to  von  Mises  yield  surface  nn+l  = - yry- 

WCn+i 

(f)  update  the  plastic  strain  £p+i  =  £p  +  A^f+1 ,  where  A£;f+1  =  A /n+lHn+l . 
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Box  1.2.  Subroutine  for  updating  damage  eigenstrain  £n+1 

Given  eigenstrain  £d  ,  damage  internal  variables  0)n ,  rn ,  microscopic  strain  £n+l ,  time  interval  At  and 
macroscopic  volumetric  strain  S  x . 

Find  £dn+u  a>n+l,  rn+l . 

If  Sn+1  >  0,  then  (tensile  loading): 

(a)  find  the  local  principal  strain  £.  by  £ 

v  '  r  r  y(«+l)  lJ(n+ 1) 

~  f  x  x>0 

(b)  compute  the  equivalent  strain  £n+l  —  \  /  \\_£,  1]  ,  where  [[xl]  =  \  ,  C 

y„+D  [ex  x<q 

is  material  parameter  indicating  damage  contribution  from  compression. 

(c)  if  e„+1  <  rn  then  r+|  =  /; ,  con+l  =  (On 
else 


V  = 

rn+\ 


A  t  . 

rn  +-^£n+\ 

& 


>  ®„+l  =®»+(r»+1-rJ 


3O(i„+0 


(d)  calculate  the  damage  eigenstrain  £n+x  =  ®„+1£„+1 . 
Else  (compressive  loading): 


Box  1.3.  Subroutine  for  updating  eigenstrain  £cn+1 

Given  eigenstrain  £c ,  lock-in  internal  variables  Sn ,  microscopic  strain  £  x ,  and  macroscopic  volumetric 


Find  £cn+l,  sn+l. 

If  Sn+1  >  0,  then  (tensile  loading): 

Sn+l=Sn’  Cl=<- 

Else  (compressive  loading): 


find  the  local  principal  strain  £,  by  £ 

1(n+ 1)  J  lJ(n+l) 

compute  the  characteristic  strain  £  n+ 1  ^ /=1  ^  £ 

if  £n+l<Sn  then 

^n+l  =  Sn  ’  £n+l  =  Sn 


^n+l  £n+l’  £n+\  ^>(.£n+ 1) 


12 


4.  Numerical  Example:  Penetration  on  Concrete  Target 

4.1.  Unit  cell  of  concrete 


penetration  model  unit  cell  model  (aggregate) 


Figure  6.  Construction  of  the  unit  cell  model  from  digital  image  of  concrete  sample 

We  start  with  a  construction  of  the  unit  cell  model  of  concrete.  The  concrete  type  is  given  in 
Section  2.  We  assume  that  crushed  limestone  constitute  the  aggregate  phase  whereas  the 
remaining  material  composes  the  mortar  phase.  The  volume  fraction  of  aggregate  is  42%  and  the 
largest  aggregate  size  is  approximately  10mm.  The  geometry  of  the  unit  cell  is  reconstructed 
from  scanned  images  in  various  cross-sections  of  the  concrete  sample.  The  finite  element  mesh  is 
then  generated  from  the  grey-level  digital  images  [30].  Figure  6  shows  the  unit  cell  model  used 
in  the  present  study.  A  thin  mortar  phase  is  formed  to  encompass  the  unit  cell  so  that  periodic 
boundary  conditions  could  be  imposed  on  the  mortar  phase  only  as  shown  in  Figure  7. 
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Figure  7.  Finite  element  mesh  of  unit  cell  model  for  concrete 


4.2.  Prediction  of  projectile  penetration 


The  numerical  model  was  implemented  in  ABAQUS  [31].  The  constitutive  law  of  concrete 
was  coded  using  Fortran  subroutines  in  UMAT/VUMAT  of  ABAQUS.  Material  parameters  were 
calibrated  using  laboratory  tests  [8].  We  adopted  the  test  data  set  of  concrete  from  [32]  for  the 
purpose  of  comparison.  The  length  and  diameter  of  the  specimens  were  on  average  111.37mm 
and  50.4mm,  respectively.  We  used  uncon  lined  compression  and  uniaxial  compression  tests  to 
calibrate  material  parameters  for  each  phase.  The  responses  for  the  two  tests  are  shown  in  Figure 
8  and  Figure  9.  The  calibrated  microscopic  material  parameters  are  listed  as  in  Table  3. 

Table  3.  Calibrated  microscopic  material  parameters  of  concrete 


Phase 

E  (MPa) 

V 

(MPa)  h  (MPa) 

e 

Mortar 

40000.0 

0.2 

11.0 

2000 

1 

Aggregate 

60000.0 

0.2 

75.0 

4000 

1 

Phase 

S  (MPa) 

G  ( MJ t m2) 

c 

P, 

Pf  ^ max 

3 

Mortar 

1.0 

1.25  xl0~5 

0.0 

0.0025 

0.1450.052 

0.05 
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Figure  8.  The  response  comparison  in  unconfined  compression 


Figure  9.  The  response  comparison  of  uniaxial  compression 


With  material  parameters  for  each  phase  defined,  we  constructed  the  finite  element  model 
for  concrete  target  and  projectile.  The  finite  element  size  in  the  penetrating  region  was  10mm. 
The  projectile  was  fired  at  five  different  velocities,  which  are  238.1 mis,  215.1ml s, 
314/« /s,  369.5/72/ 5,  and  456.4/77  / s  .  Two  types  of  properties  were  investigated,  penetration 
depth  and  projectile  deceleration  [32].  The  depth  comparison  between  the  test  data  and 
simulation  are  listed  in  Table  4.  Figures  10-14  show  the  deceleration  comparison  for  different 
striking  velocities  Vs . 

Table  4.  Penetration  depth  comparison  for  1.83m  diameter  concrete  targets 


Striking  Velocity 

(m/s) 

Penetration  Depth 

(m) 

Simulated  Depth 

(m) 

238.1 

0.30 

0.284 

275.7 

0.38 

0.365 

314.0 

0.45 

0.457 

369.5 

0.53 

0.60 

456.4 

0.94 

0.89 
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[X1.E6] 


Time  (s) 


Figure  10.  Projectile  deceleration  comparison  under  V  =  238. 1 m  /  s 


Figure  11.  Projectile  deceleration  comparison  under  V  =  275.7 m  /  s 
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Deceleration  (m/s2)  ■>?  Deceleration  (m/s2) 


[X1.E6] 


Figure  12.  Projectile  deceleration  comparison  under  V  =  3 14m  /  s 


Figure  13.  Projectile  deceleration  comparison  under  V  =  369.5 m  /  s 
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Ixl  E6] 


Figure  14.  Projectile  deceleration  comparison  under  Vs  =  456.4m  /  s 

The  depth  and  deceleration  prediction  agreed  reasonably  well  with  the  test  data.  In  this 
study,  we  also  investigated  the  damage  parameter  in  the  mortar  phase  of  the  concrete  target  at 
various  projectile  velocities  as  shown  in  Figure  15.  It  can  be  seen  that  the  damage  level  increases 
with  increase  in  striking  velocity.  The  damage  contours  depicted  in  Figure  15  provide  the 
qualitative  shape  of  the  crater  during  the  penetration. 
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Velocity=238.1m/s 


Velocity=275.7m/s 


SDV17 
(Avg:  75%) 

+2.772e-C-5 
+2.541e-C5 
+2.310e-05 
+2.07e*-C5 
+1.&4fie-C5 
+1.617e-05 
+1.386e-05 
+1.155e-05 
+9.240e-06 
+3.93C*-C6 
+4.S2Ce-Co 
+2.310e-06 
-K).OOCe+CO 


Velocity=314m/s 


Velocity=369.5m/s 


Velocity=456.4m/s 


Figure  15.  Equivalent  damage  strain  in  the  mortar  phase  of  concrete  target  penetrated  by  projectiles  at 

various  striking  velocities 


5.  Conclusions  and  future  research  directions 


A  multiscale  approach  for  predicting  penetration  of  projectiles  into  concrete  targets  has 
been  explored  as  an  alternative  to  conventional  methods  based  on  either  empirical-analytical 
approaches,  experimental  techniques  or  numerical  approaches  employing  phenomenological  or 
rules-of-mixture  constitutive  models  of  concrete.  The  multiscale  approach  presented  in  this  paper 
is  based  on  eigendeformation-based  homogenization  which  approximates  the  solution  fields  in 
the  unit  cell  in  terms  of  the  residual-free  functions  and  consequently  avoids  costly  equilibrium 
calculations.  The  inelastic  behavior  of  concrete  is  modeled  using  pore  compaction  (or  lock-in) 
and  rate-dependent  damage  eigenstrains  in  the  mortar  phase  as  well  as  plastic  eigenstrains  in 
both  the  aggregates  and  mortar.  The  method  has  been  validated  on  a  rigid  projectile  penetration 
problem  for  which  experimental  data  exists. 

Future  studies  will  focus  on  the  following  three  issues:  (i)  consideration  of  steel 
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reinforcement,  (ii)  accounting  for  uncertainty  in  the  concrete  microstructure,  (iii)  development  of 
the  mechanistic  approach  to  pore  compaction  which  so  far  has  been  accounted  for 
phenomenologically.  To  account  for  steel  reinforcement  we  will  employ  a  three-scale  reduced 
homogenization  approach  [23,24,26]  where  the  concrete  microstructure  will  be  considered  at  the 
microscale  and  a  rebar  embedded  in  a  homogenized  concrete  will  comprise  the  mesoscale.  To 
account  for  the  uncertainty  in  the  concrete  microstructure  we  will  employ  a  combination  of  the 
reduced-order  and  stochastic  homogenization  methods.  To  account  for  pore  compaction, 
development  of  large  deformation  eigendeformation  based  homogenization  will  be  pursued. 
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Appendix  A.  Brief  Theory  of  Reduced-Order  Homogenization 


A.l.  Mathematical  homogenization 


In  the  mathematical  homogenization  theory,  various  fields  are  assumed  to  be  a  function  of 
the  macroscopic  coordinates,  x ,  and  as  the  microscopic  coordinate  system,  y  =  xl  g ,  where 
£  is  a  small  positive  parameter.  Combining  the  governing  equations  (1-3)  with  the  asymptotic 
expansions 

w,  (x,  y,t)  =  uf  (x,  t)  +  C  u)  (x,  y ,  t)  +  •  •  • 


£y  (x,  y,  t )  =  s°  A  0 + C£l  (x,  y,t)+— 

(X,  v,  t )  =  cr°  (x,  y,  t )  +  £cr‘  (x,  y,t)  +  — 
1  Mij  (x,  v,  t)  =  (x,  y,  t )  +  ?n\  (x,  y,t)-\ 


yields  the  unit  cell  problem 


Lmi(y) 


Ski' 


(x,  t )  +  uk  },  (x,  y,  0  -  X  1  Ah  (x,  y,  t ) 


=  0 


>yj 


from  which  the  macroscale  stress  can  be  computed 


-  i  f  _  i  7  ) 

(Tijjix,  0  =  —  j0  Lm  O)  sk,(x,  t )  +  u\yi  (x,  y,  0  -  X  Vh  (x,  y,  t)  dy 


(Al) 


(A2) 


(A3) 
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where  0  is  the  unit  cell  domain,  while  aj  and  £ij  are  macroscale  stress  and  strain, 
respectively. 

A.  2.  Reduced  order  multiscale  system 


Following  [23],  the  microscale  displacement  field  u)  (x,y,t)  is  constructed  so  that  the 
stress  field  in  the  unit  cell  would  automatically  satisfy  equilibrium  equations.  This  is 
accomplished  by  introducing  eigenstrain  1  jJl  and  eigenseparations  Sn 

u)(x,y,t)  =  <W,0  +  Z  + 

I 

=  0 + E  le  ‘rT  ( y<  A1  A  y,  <  A  (A4) 


+lKr(y,mXx,y,l)dy 

The  resulting  microscale  displacement  gradients  are  given  by 

u),yj  (*>  y,  0  =  Gijkl  £kl(x,  t)  +  £  j0 1  gfd  ( y,  y)'  mI,  O,  y,  t)dy 

+ ■  Js  gth  (y>  y)sn  y, 


(A5) 


where  Gjkl ,  1  ,  and  gfd’  are  influence  functions  for  macroscale  strain,  eigenstrain,  and 


eigendisplacement,  respectively,  which  can  be  computed  by  solving  a  sequence  of  elastic 
boundary  value  problems  prior  to  nonlinear  macro  analysis.  A  reduced  order  model  is  obtained 
by  discretizing  eigenstrain  and  eigenseparation  fields  as 


’  4(x,y,t)^INu\y)Ifig\x,t) 

a= 1 
m 

5h(xJd)  = 

f=i 


(A6) 


where  and  m  are  the  numbers  of  partitions  of  phases  and  interfaces  respectively,  Ijuijia) 
and  8-f  ]  are  the  average  eigenstrain  and  eigendisplacement  in  the  phase  partition  a  and 
interface  partition  £ ,  respectively.  Nia)(y)  is  a  piecewise  constant  shape  function  defined  as 


Nia\y ) 


jl  ye0(a\ 
jo  y  e  0(“\ 


(A7) 


A(f)(y) 


is  a  linear  combination  of  piecewise  linear  finite  element  shape  functions  defined  over 


partition  g . 

Combining  the  unit  cell  equations  (A2)  with  the  decomposition  (A5)  with  the  above 
discretization  yields  the  reduced  system  of  equations: 
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where  t-  =  G(gH)  represents  the  traction  along  interface.  All  coefficient  tensors  in  the 

reduced-order  system,  such  as  1  Pijkl,  Qjh ,  Aijkl,  ‘C-j,  Dnm ,  Bni] ,  and  ^... »  are 

determined  prior  to  the  nonlinear  analysis  in  the  preprocessing  stage.  The  detailed  formulation 
can  be  found  in  [23,24]. 
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